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Afe&EXISl 

Two  topics  on  the  numerical  solution  of 
boundary- integral  equations  arising  in  linear 
wave-body  interactions  are  discussed.  The 
properties  of  a  spectral  technique  for  the 
solution  of  the  integral  equation  are  analysed 
and  eosipared  to  the  conventional  collocation 
method.  It  is  shown  that,  using  this  tech¬ 
nique,  hydrodynamic  forces  predicted  by  the 
source-distribution  method  are  identical  to 
those  obtained  from  the  direct  solution  for 
the  velocity  potential.  The  second  part  of  the 
paper  investigates  the  numerical  properties  of 
s  method  which  removes  the  effects  of  the 
irregular  frequencies  for  bodies  of  general 
shape  at  a  small  cosrputational  and  algorithmic 
overhead.  Its  performance  is  illustrated  in 
the  evaluation  of  the  heave  and  sway  hydro- 
dynamic  coeff icientg_of  a  circle  and  a  rectsn- 
81'-  ■  •'  ,,  V  J  .  /  , 

1.  Introduction  > 


The  solution  of  boundary-integral  equa¬ 
tions  for  the  evaluation  of  the  linear  wave 
loads  on  marine  structures  is  a  common  task  in 
today's  practice.  Its  success  is  due  to  its 
algorithmic  simplicity,  the  ease  of  describing 
the  surface  of  a  three-dimensional  body  by  a 
collection  of  facets  and  the  moderate  size  of 
the  linear  systems  to  be  solved.  These  are 
illustrated  by  its  widespread  use  by  sero- 
dynamicists  [  Hess  and  Smith  (1966)  ]. 


Q  In  the  presence  of  a  free  surface,  the 
vJ  efficiency  of  the  method  reli«s  on  the  fast 
evaluation  of  the  wave-source  potential  which 
L.L.J  is  a  substantially  more  complex  function  to 
— I  compute  than  its  counterparts  in  an  infinite 
fluid  and  an  acoustic  medium.  Existing  methods 
for  the  computation  of  its  values  and  deriva¬ 
tives  are  hard  to  evaluate,  since  it  is  the 
performance  of  the  integrated  radiation-dif¬ 
fraction  computer  programs  that  is  usually 
■"■^5-eported.  For  the  three-dimensional  comput¬ 
ation;  reported  in  the  present  paper,  a  set  of 
very  efficient  algorithms  developed  by  Newman 
fl9F5a)  for  water  of  finite  and  infinite  depth 
and  coded  in  the  subroutine  FINGREEN  have  been 
utilized. 


cjsL 


A  distinct  feature  of  wave  boundary- 
integral  equations  are  the  "  irregular  " 
frequencies.  They  coincide  with  the  eigen- 
frequencies  of  the  interior  Dirichlet  or 
Neumann  problems ,  and  are  known  to  introduce 
large  errors  in  the  predicted  hydrodynamic 
forces,  often  over  a  substantial  band  of 
frequencies.  A  comprehensive  analysis  of  the 
mathematical  properties  of  boundary- integral 
equations,  (with  emphasis  in  acoustics),  along 
with  a  survey  of  techniques  used  for  the 
removal  of  the  irregular  frequencies,  is  given 
in  the  recent  book  of  Colton  and  Kress  (1983). 
The  numerical  aspects  of  boundary-integral,  as 
well  as  finite-element,  hybrid-integral  and 
f inite-element /boundary-integral  methods  in 
free-surface  flows  are  reviewed  by  Mei  (1978), 
Yeung  (1982)  and  Euvrard  (1983). 

The  first  part  of  the  paper  analyzes  the 
properties  of  a  technique  for  solving  bound- 
' ary-integral  equations.  It  is  often  quoted  in 
/>~the  literature  as  the  Galerkin  method.  In  most 
implementations  of  the  boundary-integral  for¬ 
mulation,  the  body  surface  is  approximated  by 
a  collection  of  H  plane  quadrilaterals  or  tri- 
sngles.  The  conventional  collocation  method 
enforces  the  validity  of  the  equation  at  a 
single  point  on  each  panel,  usually  the  cen¬ 
troid.  'hie  method  proposed  here,  averages  the 
equation  over  each  panel  and  avoids  the  need 
to  select  a  collocation  point.  It  belongs  in 
the  general  category  of  "spectral"  tech¬ 
niques  which  express  the  solution  in  terms  of 
N  basis  functions,  and  then  pre-integrate  the 
product  of  the  equation  to  be  solved  with  each 
function  of  this  set.  Here,  the  i-th  basis 
function  has  the  value  one  on  the  i-th  panel 
and  zero  on  the  rest  of  them. 

The  Galerkin  technique  has  a  set  of  inter¬ 
esting  symmetry  properties.  The  diagonal  ele- 
ements  in  the  added-mass  and  damping  matrices 
obtained  from  the  source-distribution  and  the 
Green  method  are  identical.  The  off-diagonal 
coefficients  A..,  B .  .  obtained  from  one 
method  are  ideniicai^to  the  A..,  B..  coef¬ 
ficients  which  follow  from  ^ the  Vther.  More¬ 
over,  the  exciting  forces  obtsined  from  the 
solution  of  the  Green  integral  equation  for 
the  diffraction  velocity  potential,  are 
identical  to  those  obtained  from  the  Haskind 
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relation*  with  the  radiation  velocity  poten-  degree*  out  of  pbaae  from  the  velocity  poten¬ 
tial  supplied  by  the  aour ca-diatribut ion  tial,  or  in  phaae  with  the  preaaure.  Thua 

aethod,  and  vice  veraa.  Analogou*  reaulta  do  energy  flow  out  of  the  interior  doaain 

not  hold  in  the  collocation  aethod.  preventing  the  peraiatence  of  eigenaolutiona . 

The  direction  of  the  energy  flux  ia  controlled 

Coaputationa  of  the  hydrodynaaic  coef-  by  the  aign  of  the  conatant  uaed  in  the  linear 

ficient*  for  a  aeai-aubaerged  apheroid  and  a  coabination  of  the  two  equation*, 

truncated  vertical  circular  cylinder  have  been 

carried  out  by  Breit,  Newman  and  Sciavouno*  Euvrard,  Jami,  Lenoir  and  Martin  (1981) 

(1985).  No  aubatantial  difference  from  the  were  the  firat  to  adapt  thia  methodology  to 

prediction*  of  the  collocation  aethod  ha*  been  wave-body  interaction  problem*,  their 

obaerved.  Near  the  irregular  frequenciea,  the  formulation  coabinea  a  layer  of  finite 

Galerkin  aethod  led  to  a  reduction  of  the  element*  in  the  fluid  doaain  which  encloaea 

error  and  the  frequency  bandwidth  over  which  the  body  boundary,  with  a  boundary- integral 

it  occur*.  In  principle,  it  require*  an  repreaentation  in  the  exterior  doaain  analo- 

additional  numerical  integration  for  the  gous  to  that  outlined  in  the  preceding  para- 

evaluation  of  the  infuence  of  panel  i  on  panel  graph.  Computation*  of  the  hydrodynaaic  coef- 

j.  In  Section  2,  an  algorithm  ia  propoaed  ficienta  of  tbree-diaenaional  bodie*  were 

which  reduce*  aubatantially  thia  overhead,  found  free  of  the  effect*  of  irregular 

while  maintaining  the  main  feature*  of  the  frequenciea. 

Galerkin  technique. 

In  the  preaent  paper  the  method  of 

For  tbe  Green  integral  equation,  the  ir-  Burton  and  Miller  (1971)  i*  applied  to  the 

regular  frequenciea  coincide  with  the  eigen-  aolution  of  the  radiation  problem.  It  corre- 

frequenciea  of  the  interior  Dirichlet  problem.  aponda  to  the  limit  of  a  finite-element  layer 

Both  in  two  and  three  dimenaiona ,  they  can  be  of  zero  thickneaa  in  the  acheae  of  Euvrard  et 

auppreaaed  by  adding  an  artificial  lid  on  the  «1*  A  circle  and  a  rectangle  have  been  anal- 

interior  free  aurface  a*  auggeated  by  Ohmatau  yzed.  Their  boundary  haa  been  approximated  by 

(1975)  in  connection  with  the  aource-diatr i-  atraight  segments  and  the  equation  has  been 

but  ion  method,  and  by  Kleimnan  (1982)  for  the  aatiafied  at  their  aid-point,  according  to  the 

Green  integral  equation.  Thia  approach  ia  collocation  aethod.  In  both  the  heave  and  away 

effective,  but  may  lead  to  a  aubatantial  added  maas  and  damping  coefficients,  errors 

increase  in  the  computational  effort,  cape-  have  been  suppressed  at  and  in  the  vicinity  of 

cially  in  three  dimenaiona.  Ogilvie  and  Shin  the  irregular  frequencies. 

(1977)  suggested  an  alternative  route  by 

adding  a  wave  source  at  the  origin  of  the  The  associated  computational  and  algo- 

coordinate  system,  assumed  to  lie  on  the  rithmic  overhead  is  amall,  aince  the  effort 

interior  free  surface,  acting  a*  an  absorber  involved  in  the  computation  of  the  second 

of  the  energy  associated  with  the  interior  derivatives  of  the  wave  source  potential  (they 

Dirichlet  eigenaolutiona.  Thia  approach  was  appear  in  the  kernel  of  the  equation  of  the 

implemented  in  two-dimenaiona  and  was  shown  to  first  kind)  is  not  large  relative  to  that 

successfully  remove  the  firat  irregular  required  by  its  value  and  firat  derivatives, 

frequency,  at  a  small  additional  computational  In  deep  water  thia  overhead  is  negligible 

effort.  Ursell  (1982)  later  established  that  because  of  the  existence  of  recurrence 

any  number  of  irregular  frequencies  can  be  relation*  which  relate  derivatives  of  high  to 

removed  if  a  sufficient  number  of  singular-  those  of  lower  order.  Moreover,  the  size  of 

ities  are  added  at  the  origin.  No  implements-  the  linear  system  is  unaffected  by  the 

tion  of  this  method  in  three  dimensions  is  superposition  of  two  equations  over  the  same 

known  to  the  authors.  number  of  panels. 

Related  work  in  acoustics  predate*  the  The  aethod  is  currently  being  extended  to 

studies  of  marine  hydrodynamicists  by  about  a  three-dimensional  problems  where  the  irregu- 

decade.  Refe  rences  to  early  atudiea  can  be  lar-f requency  effects  are  generally  less 

found  in  Colton  and  Kress  (1983).  Burton  and  pronounced,  and  is  expected  to  be  effective 

Miller  (1971)  proposed  a  methodology  which  for  bodies  of  general  shape.  This  ia  believed 

exploits  the  different  location  of  the  irregu-  to  b*  the  case  because  the  addition  of  the 

lar  frequencies  of  integral  equations  of  the  integral  equation  of  the  first  kind  to  the 

first  and  second  kind.  The  linear  combination  Green  equation  essentially  corresponds  to  an 
of  two  such  equations  for  the  exterior  Neumann  additional  singularity  distribution  on  the 
problem  haa  a  unique  solution  on  the  entire  actual  body  surface  rather  than  at  a  selected 

real  frequency  axis,  and  thus  ia  free  of  ir-  »*t  of  points  internal  to  its  boundary, 

regular  frequencies.  The  condition  for  thia  to 
occur  is  that  one  of  the  two  equations  must  be 
multiplied  by  the  imaginary  unit  times  a 
positive  conatant.  It  turns  out  that  the  asso¬ 
ciated  interior  homogeneous  problem  ia  of 
mixed  Dirichlet  and  Neumann  type,  the  two 
being  90  degree*  out  of  phase.  Non-trivial 
eigen  solutions  cannot  exist  since  on  the 
interior  boundary  the  normal  velocity  ia  90 


We  ere  interested  to  evaluate  the  hydro- 
dynamic  pressure  on  the  surface  of  a  body 
interacting  with  regular  free-surf ace  waves. 
Linearity,  and  the  assumption  of  irrotational 
flow,  allows  the  reduction  of  the  problem  to 
the  evaluation  of  a  velocity  potential  $  (z) 
which  satisfies  the  Laplace  equation  in  the 
fluid  domain,  the  linear  free-surface  con¬ 
dition,  a  prescribed  normal  gradient  v(x)  on 
the  body  surface,  the  vanishing  of  its  gradi¬ 
ent  at  large  depths  and  a  radiation  condition 
at  infinity. 

Tvo  boundary-integral  formulations  are 
common,  and  both  can  be  derived  from  Green's 
theorem.  According  to  the  source-distribution 
method,  $(x)  is  represented  by  a  distribution 
of  wave  sources  on  the  body  surface, 


*(x)  -  /  a(£)Cfx;5)de  (2.1) 

S 

where  G(x;  5  )  is  the  wave-source  potential  at 
the  point  x  due  a  unit  source  located  at  the 
point  i  .  Enforcing  the  normal  velocity  on 
the  body  boundary  leads  to  an  integral  equa¬ 
tion  for  the  strength  5  (x)  of  the  source 
distribution, 


angles,  assuming  that  the  unknown  functions 
have  constant  values  on  each  panel.  The 
equations  are  enforced  at  a  set  of  collocation 
points  usually  taken  to  be  the  centroids.  The 
Galerkin  technique,  proposed  here,  averages 
instead  the  equations  over  esch  panel.  In 
discrete  form,  equations  (2.1)-(2.2)  become. 


1  N 

0.  -  7-  I  «  J  dx  /  d£  G(x;f 
Ai  j-1  J  S .  S, 


Ai  j-1 


G  -  J  dx  J  d£  G(L;x) 
iJ  S.  S, 
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and  equation  (2.3), 
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The  application  of  a  different  variant  of 
Green's  theorem  leads  to  an  integral  equation 
with  the  velocity  potential  itself  as  the 
unknown , 


1  ..  *  3G(ti  ,x) 

-  -  j(x)  +  ;  iU)  — ^ —  d£ 


:  G^;x)dc 


Equations  (2.2)  and  (2.3)  are  adjoint  Fredholm 
integral  equations  of  the  second  kind,  since 
the  kernel  of  (2.2)  is  obtained  from  the 
transposition  of  the  arguments  of  the  kernel 
of  (2.3).  The  preference  of  the  one  versus  the 
ocher  depends  on  the  application  for  which 
they  are  being  considered.  If,  for  example, 
the  flow  velocities  sre  required  then  equs- 
tions  (2.1 )— ( 2.2)  sppesr  to  be  at  an  advantage 
since  the  evslustion  of  the  second  derivstives 
of  G  is  not  necessary.  If,  on  the  other  hand, 
only  quantities  dependent  on  the  velocity  po- 
tentisl  are  needed,  then  (2.3)  must  be  prefer¬ 
red  due  to  the  reduced  storage  requirements. 

Their  numerical  solution  is  usually  ob¬ 
tained  by  approximating  the  body  boundary  by  a 
collection  of  plane  quadrilaterals  or  tri- 


J  dx  J  dS 


where  A.  is  the  area  of  the  i-th  panel.  The 
integration  with  respect  to  the  x-variable 
introduced  by^^ie  Gherkin  averaging,  allows 
the  matrices  D  ,  D  and  G  to  preserve 
the  symmetry  properties  of  the  corresponding 
operators  in  the  continuous  case.  In  partic¬ 
ular,  the  matrix  G  is  symmetric,  thus 


Gij  '  Gji 


(2.10) 


and  the  matrices  and  sre  the 


transpose  of  each  other,  or 


D(G)  -  D(S> 
ij  Ji 


(2.11) 


The  proof  of  (2.10)  follows  from  the  symmetry 
of  the  wave-source  potential  with  respect  to 
its  arguments,  and  of  (2.11)  by  exchanging  the 
role  of  the  dunssy  variables  under  the  integral 
signs  in  (2.7)  and  transposing  the  i  and  j 
indices.  Analogous  results  do  not  hold  in  the 


collocation  method,  vh«r«  the  corresponding 
matrix  elements  are  obtained  by  replacing  the 
integration  with  respect  to  the  x-variable  in 
(2.4)-(2.9)  by  the  selection  of  a  collocation 
point. 

Let  A  *  diag(  A J ,  and  define, 

D  *  •  j  A  +  D(C)  (2.12) 

The  solution  for  the  velocity  potential 
obtained  from  equations  (2.4)-(2.7)  in  terms 
of  the  matrices  A,  D  and  G,  ia  given  by 

/  c'j  —1  T  -1  ** 

♦  W  -  A  GCD1)  A  v  (2.13) 


Both  H  and  B  are  complex  scalars. 
Three  properties  follow  from  equations  (2.16) 
and  (2.17)  : 


1)  For  the  diagonal  hydrodynamic  coefficients 

Sk"  i  Bkk/  “  .k  6,  the 

vectors  u.  and  v.  are  equal  to  the 
values  on^  the  j-^h  panel  of  the  unit  vector 
n^  which  points  out  of  the  fluid  domain. 
In  this  case  the  hydrodynamic  coefficients 
predicted  by  the  two  methods  are  identical, 
since 


7  S  v  ■  7  #T  v  (2.18) 


and  the  solution  of  (2.8)  by 


2)  For  the  off-diagonal  coefficients,  it 
follows  from  (2.16)  and  (2.17)  that 


;(G)  -  D_1G  7  (2.14) 


For  an  arbitrary  normal-velocity  vector  v,  a 
necessary  condition  for  the  identity  of  the 
two  solutions  (2.13)_^nd  (2.14)  is  the  equa¬ 
lity  of_thj  matrix  D  G  with  the  matrix 
A  G  (DT)  A,  or  equivalently  the 
symmetry  of  the  matrix 


■  »-(C> 


IK 


(2.19) 


3)  For  the  diffraction  exciting  force,  we 
define 
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■(  — -  ) 
v  3n  'j 


(2.20) 


W  -  A  D_1G  (2.15) 


A  proof  that  H  is  symmetric  did  not  prove 
possible.  Numerical  experiments  for  a  model 
problem  in  two  dimensions  in  an  infinite  fluid 
revealed  that  the  matrix  H  is  "almost  synme- 
trie",  meaning  that  elements  with  symmetric 
locations  relative  to  the  principal  diagonal 
agreed  to  2-3  significant  digits.  This  sug¬ 
gests  the  proximity  of  the  values  for  the 
velocity  potential  obtained  from  each  method. 

The  hydrodynamic  forces  can  be  obtained 
from  the  solutions  (2.13)  and  (2.14)  by  multi¬ 
plying  the  velocity  potential  by  the  panel 
area  A.  and  the  vector  u.  which 
represents  the  "direction"  of  the  force  we  are 
interested  to  evaluate.  This  operation  is 
equivalent  to  a  pre-multiplication  of  tbe 
velocity-potential  vector  by  the  vector 
(Au)  .  The  resulting  hydrodynamic  force 
obtained  form  the  source-distribution  method 
is  given  by 


HlS'  »  u‘w‘  v  (2.16) 


and  from  the  direct  solution  of  the  Green 
integral  equation, 


where  ? .  ii  the  incident-wave  velocity 
potential.  IfX^is  the  diffraction 
exciting  force  in  the  k-th  direction,  it 
can  be  deduced  from  (2.16)  and  (2.17)  that 
the  force  predicted  by  the  source-distri¬ 
bution/Green  method  by  directly  solving  the 
diffraction  problem,  is  identical  to  the 
force  obtained  from  the  use  of  the  Haskind 
relation  with  the  velocity  potential 
supplied  by  the  solution  of  the  Green/ 
source-distribution  integral  equation. 

Computations  of  the  hydrodynamic  coeffi¬ 
cients  of  a  spheroid  and  a  vertical  circular 
cylinder  using  both  methods  have  been  carried 
out  by  Breit,  Newman  and  Sclavounos  (1985).  A 
radiation-diffraction  computer  program  has 
been  written  for  *.he  hydrodynamic  analysis  of 
bodies  of  general  shape.  Their  wetted  surface 
is  approximated  by  a  collection  of  plane 
quadrilaterals  and  triangles,  as  illustrated 
in  Figures  1  and  2  for  a  quarter  of  the  spher- 
roid  and  the  vertical  cylinder  respectively. 

For  inter-panel  distances  small  compared 
to  their  characteristic  dimensions,  the 
Rsnkine  singularity  (including  when  appropri¬ 
ate  its  image  with  respect  to  the  free  surface 
and  the  sea  bottom)  is  subtracted  from  the 
wave-source  potential  and  integrated  analytic¬ 
ally  over  the  panels.  For  large  distances 
betveen  the  panels  the  total  wave-source 
potential  is  integrated  by  quadrature.  The 
algorithms  for  the  integration  of  the  Rankine 
singularities  on  plane  quadrilaterals  and  for 
the  evaluation  of  the  wave-source  potential 


(2.17) 


have  bean  developed  by  Neman  (1985a  &  b)  and 
coded  in  Che  subroutines  FINGREEN  and  RPAN 
respectively .  The  four-node  Cauas-Legendr e 
quadrature,  adapted  to  a  plana  quadrilateral 
of  general  shape,  has  been  used  for  the 
evaluation  of  influence  coefficients  in  both 
the  collocation  and  Galerkin  aethoda.  Sugges¬ 
tions  on  the  order  of  the  quadrature  to  ibe 
used  in  a  production  radiation-diffraction 
coaputer  program  are  discussed  at  the  end  of 
the  section. 

The  hydrodynaaic-force  predictions  of  the 
collocation  and  Galerkin  methods  did  not 
differ  substantially  away  froa  the  irregular 
frequencies.  Tabulated  results  of  high 
accuracy  are  reported  in  Breit,  Neman  and 
Sclavounos.  Figures  3  and  4  illustrate  the 
behaviour  of  the  two  methods  at  the  first 
heave  irregular  frequency  of  the  spheroid  and 
the  cylinder  respectively.  The  solid  lines 
represent  the  predictions  of  an  independent 
curvilinear-panel  program  for  the  spheroid, 
and  a  Fourier-transformed  time-domain  solution 
for  the  cylinder.  The  Galerkin  predictions 
appear  to  be  less  sensitive  to  the  irregular- 
frequency  errors,  especially  for  the  spheroid 
coefficients. 

The  Calerkin  technique  requires  an  addi¬ 
tional  integration  for  the  evaluation  of  each 
element  in  the  D.  .  and  G..  matrices,  [eq. 
(2.4),  (2.9))  Vtlative^  to  the  collocation 
method.  It  is  reasonable  to  assume  that  the 
accuracy  in  the  integration  over  the  i-th 
panel  needs  to  be  no  higher  than  that  over  the 
j-th  panel.  Concerning  the  Rankine  source  and 
dipole,  analytical  expressions  for  the  double 
integral  over  a  pair  of  plane  quadrilaterals 
are  not  known  to  the  authors.  When  the  analyt¬ 
ical  expressions  are  utilized  for  the  evalu¬ 
ation  of  the  Rankine  source  and  dipole  inte¬ 
grals  over  the  j-th  panel,  a  four-node  Gauss- 
Legendre  quadrature  is  suggested  for  the  inte¬ 
gration  over  the  i-th  panel.  Since  this  result 
is  frequency- independent  it  may  be  evaluated 
once  and  stored. 

The  integrals  of  the  remaining  slowly- 
varying  but  frequency-dependent  parts,  can  be 
evaluated  using  a  quadrature  scheme  of  the 
same  order  for  the  j-th  and  i-th  panels.  In 
the  collocation  method,  the  use  of  a  four-node 
Gauss  quadrature  causes  an  increase  by  a 
factor  of  four  in  the  number  of  evaluations  of 
the  wave-source  potential,  versus  the  single¬ 
node  centroid  integration.  This  factor  may  be 
as  high  as  sixteen  in  the  Galerkin  method. 
Optimality  requires  that  errors  due  to  qua¬ 
drature  and  the  approximation  of  the  geometry 
by  plane  panels  must  be  of  comparable  mag¬ 
nitude.  This  may  be  achieved  by  increasing  the 
number  of  panels  and  utilizing  a  single-node 
quadrature.  This  decision  depends  on  the 
efficiency  in  the  evaluation  of  the  wave- 
source  potential,  the  solution  of  the  linear 
system  and  the  computing  environment.  If  the 
single-node-quadrature  route  is  selected,  the 
collocation  and  Galerkin  methods  are  compara¬ 
bly  expensive  over  s  large  number  of  fre- 
quenc ies . 


Figure  1  :  Discretization  of  a  quarter  of  the 

surface  of  a  prolate  sphepoid  (B/L  ■ 
1/8)  by  64  panels. 


Figure  2  :  Discretization  of  s  quarter  of  the 
wetted  surface  of  a  truncated 
vertical  cylinder  (R/T  ■  2)  by  128 
panels . 


Figure  3  :  Heave  added-mass  and  damping 
coefficients  of  a  prolate  spheroid  (B/L  -  1/8) 
near  the  first  irregular  frequency,  made  non- 
dimensional  by  the  diaplaced  volume  the  water 
density  and  the  frequency  of  oscillation.  The 
solid  line  is  obtained  from  an  independent 
curvilinear-panel  program  and  the  tick  marks 
are  predictions  form  the  plane-panel  program. 
[From  Breit,  Newman  and  Sclavounos  (1983)] 


Figure  4  :  Heave  added-maas  and  damping 
coefficients  of  s  truncated  vertical  cylinder 
near  the  first  irregular  frequency,  made  non- 
dimensional  by  the  displaced  volume  the  water 
density  and  the  frequency  of  oscillation.  The 
solid  line  is  obtained  from  the  Fourier  trans¬ 
form  of  an  independent  time-domain  program  and 
the  tick  marks  are  predictions  form  the  plane- 
panel  program.  [From  Breit,  Nevman  and 
Sclavounos  (1983)] 


The  Green  integral  equation  (2.3)  is  known 
to  possess  non-vanishing  homogeneous  solutions 
at  a  discrete  set  of  frequencies  which  corre¬ 
spond  to  the  eigenf requeue ies  of  the  interior 
Dirichlet  problem.  Their  detrimental  effect  in 
the  numerical  predictions  of  the  added-masa 
and  damping  coefficients  of  surface-piercing 
bodies  has  been  illustrated  in  Figures  3  and 
4. 

Although  a  discrete  set  in  the  continuous 
ease,  their  presence  in  the  discrete  problem 
is  manifested  by  substantial  errors,  often 
over  a  quite  vide  frequency  band  around  their 
exact  location.  This  is  due  to  the  bad  "con¬ 
ditioning"  of  the  integral  equation  (2.3)  not 
only  at,  but  also  in  the  vicinity  of  the  ir¬ 
regular  frequencies.  Bad  conditioning  is 
known  to  cause  large  errors  in  the  solution 
when  a  small  perturbation  is  imposed  on  the 
equation.  In  wave-body  interactions  sources  of 
such  perturbations  are: 

1)  The  approximation  of  the  body  boundary. 

2)  The  approximation  of  the  velocity  poten¬ 
tial  . 

3)  Errors  in  the  evaluation  of  the  wave- 
source  potential. 

4)  Quadrature  errors  in  the  evaluation  of 
the  influence  coefficients. 

5)  The  approximate  way  in  which  the  equation 
is  being  satisfied. 

6)  Roundoff  errors  in  the  solution  of  the 
linear  system. 


A  measure  of  the  ratio  of  the  output 
versus  the  input  errors  in  the  solution  of 
integral  equations  is  often  supplied  by  the 
"condition  number".  Explicit  definitions  of  it 
are  known  for  matrix  equations.  Thus  the 
discrete  form  of  an  integral  equation  may  be 
used  to  obtain  an  estimate  of  it. 

Numerical  experiments  indicate  that  the 
errors  and  frequency  bandwidth  of  the  irregu¬ 
lar  frequencies  decrease  with  increasing 
numbers  of  panels.  The  associated  computation¬ 
al  cost,  however,  prevents  this  to  be  consid¬ 
ered  s  viable  treatment  in  practice.  A  short 
survey  of  successful  methods  for  the  removal 
of  the  irregular-frequency  effects  has  been 
given  in  the  Introduction. 

The  method  of  Burton  and  Miller  (1971), 
developed  for  the  solution  of  an  acoustic 
scattering  integral  equation,  it  here  adapted 
to  the  wave-body  interaction  problem.  The 
Green  equation  (2.3)  is  valid  for  a  point  x  on 
the  body  boundary.  If  x  lies  in  the  fluid 


domain,  the  factor  1/2  which  multiplies  the 
firat  term  needa  to  be  replaced  by  unity. 
Taking  the  derivative  of  both  aidea  in  the 
direction  of  the  unit  vector  n  which  point  a 
out  of  the  fluid  domain,  and  letting  the  point 
x  approach  the  body  boundary,  we  obtain 


The  numerical  conditioning  of  equation 

(3.1)  ia  worse  than  that  of  the  Green  equation 
(2.3).  Thie  ia  generally  known  to  be  true  for 
integral  equationa  of  the  firat  kind.  Bence, 
the  perturbations  1-6  cause  errors  in  the 
solution  of  (3.1)  large  relative  to  those  in 
the  solution  of  (2.3). 
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For  a  prescribed  normal  velocity,  (3.1)  ia  an 
integral  equation  of  the  firat  kind  for  the 
velocity  potential  on  the  body  boundary.  Its 
irregular  frequencies  correspond  to  the  eigen- 
frequencies  of  the  interior  Neumann  problem. 
Burton  and  Miller  (1971)  show  that  the  linear 
combination  (2.3)  ♦  i  2(3.1),  or 
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has  no  real  irregular  frequencies  for  real  and 
positive  values  of  the  parameter  2  .  numeri¬ 

cal  experimentation  suggests  that  for  *ili  . 
of  the  parameter  3  ranging  from  0.2  to  0.3, 
the  performance  of  (3.2)  is  optimal.  A 
discussion  of  the  properties  of  aquations 

(2.3),  (3.1)  and  (3.2)  for  acoustic-wave  prob¬ 
lems,  is  given  by  Colton  and  Kress  (1983). 

Approximating  the  body  boundary  by  a 
col  lection  of  panels  S.,  and  satisfying 
(3.2)  at  collocation*  points  located  at  the 
panel  centroids,  we  obtain  the  discrete  set  of 
equat ions 

.  A.  *  +  z  >  '  (1  +  ia  - — ■—  )G(£;x  )dt 

-  j  s  ’V  1 

.1  *  J  i  1 


1 . S 


(3.3) 


The  integral  of  the  double  normal  deriva¬ 
tive  of  the  singular  part  of  the  wave-source 
potential  needs  careful  interpretation  For 
i«j,  it  is  equal  to  the  normal  velocity  on  the 
panel  due  to  a  distribution  of  dipoles  of  con¬ 
stant  strength  on  its  surface.  This  value  is 
known  to  exist  at  points  not  lying  on  its 
edges. 


The  effect  of  equation  (3.1)  on  the  solu¬ 
tion  of  (3.2)  is  controlled  by  the  magnitude 
of  the  positive  quantity  a  .  If  a  *0,  (3.2) 
reducea  to  the  Green  equation.  For  finite 
values  of  a  ,  (3.2)  is  expected  to  be  free 
of  irregular  frequency  effects.  This  turns  out 
not  to  be  the  caae  in  practice  for  very  small 
valuea  of  a  which  damp  excessively  the 
effect  of  equation  (3.1).  In  this  case  the 
error  in  the  hydrodynamic  hydrodynamic  coef¬ 
ficients  near  the  irregular  frequencies  of 

(2.3),  although  reduced,  is  still  substantial. 
For  large  values  of  a  ,  the  effects  of 
equation  (3.1)  are  dominant.  The  predicted 
coefficients  are  now  erroneous  near  its  own 
irregular  frequencies.  Smaller,  but  still 
noticeable,  errors  are  alao  present  for  all 
frequencies  due  to  its  poor  conditioning.  The 
magnitude  of  both  can  be  obviously  controlled 
by  the  selection  of  a  sufficiently  small  value 
of  a  which  strikes  a  balance  between  the 
errors  coming  from  the  irregular  frequencies 
of  the  Green  equation  (2.3)  and  those  coming 
from  the  ill-conditioning  of  equation  (3.1). 

Computations  of  the  heave  and  sway  coef¬ 
ficients  of  a  circle  and  a  rectangle  (B/T-2), 
indicate  that  the  value  of  2  -0.2  produces 
satisfactory  results  over  a  wide  range  of  fre¬ 
quencies.  A  value  of  less  than  one  is  not 
surprising.  If  the  error  in  the  solution  of 

(3.2)  resulting  from  the  presence  of  equations 

(2.3)  is  to  be  comparable  to  the  error  due  to 
equation  (3.1),  the  value  of  a  must  be 
comparaable  to  the  ratio  of  their  condition 
numbers  which  is  a  quantity  with  magnitude 
less  than  one. 

A  computer  program  has  been  written  for 
the  hydrodynamic  analysis  of  sections  of  gene¬ 
ral  shape  in  deep  water  [  Sclavounos  (1985)  ). 
The  computational  effort  involved  in  the  set¬ 
up  and  solution  of  the  discrete  equations 

(3.3) ,  is  for  all  practical  purposes  compara- 
ble  to  that  required  when  01  -0.  This  is  due 
to  the  existence  of  recurrence  relations 
which  express  higher  derivatives  of  the  wave- 
source  potential  in  terms  of  derivatives  of 
lower  order.  The  size  of  the  linear  system  is 
unaffected  by  the  addition  of  (3.1)  to  the 
Green  equation. 

Figure  5  presents  the  heave  hydrodynamic 
coefficients  of  the  circle  and  the  rectangle 
near  the  first  irregular  frequency  of  equation 

(2.3) .  Predictions  for  5  “0,  0.02  and  0.2  are 
compared  to  those  obtained  from  the  hybrid- 
integral  solution  of  Nestegard  and  Sclavounos 
(1984)  which  is  known  to  be  free  of  irregular- 
frequency  effects.  The  corresponding  results 
for  sway  are  presented  in  Figure  6  with 
analogous  conclusions. 


Rectangle  )  A22  I  Rectangle  )  B22  (  Circle  )  A22  (  Circle 


822  I  Rectangle  )  A22  <  Rectangle  )  B22  I  Circle  )  A22  I  Circle  ) 

,J> .  00  0.03  0  06  0.09  ^.40  1.45  1.50  1.55  J). 00  0.05  0.10  0.15  ^  75  O.BS  0.95 


Figure*  7  »nd  8  illu*tr*te  the  performance  Euvrard,  D.,  Jami,  A.,  Lenoir,  M.  and  Martin, 

of  equation  (3.3)  near  the  fir*t  irregular  D.  (1981).  "Recent  Progre**  Toward*  an  Optimal 

frequency  of  equation  (3.1).  Re*ult»  are  pre-  Coupling  of  Finite  Element*  and  Singularity 

•ented  for  c.  «  0,  0.2  and  from  the  solution  Distribution  Procedure*  in  Numerical  Ship 

of  equation  (3.1)  alone.  The  prediction*  for  Hydrodynamic*".  Proc.  3rd  Inti.  Conf.  Num. 

oi *0  are  accurate  over  that  frequency  range,  Ship  Hydrodyn.,  Paris, 

and  are  in  good  agreement  with  the  results 

for  oi  -0.2.  The  larger  discrepancies  occur  Euvrard,  D.  (1983).  "Sur  Quelquea  Kethodes  de 

in  the  heave  damping  coefficient*,  but  these  Resolution  De*  Problemes  Lineaires  en  Dotnaine 

are  probably  due  to  its  small  values.  The  non  Borne".  ENSTA  publication, 

predictions  from  equation  (3.1)  display  the 

escpected  error  near  its  first  irregular  fre-  Hess,  J.  L.  and  Smith,  A.  M.  0.  (1966). 

quency.  Evident  is  also  a  non-negligible  er-  "Calculation  of  Potential  Flows  Around 

tot  over  a  wider  frequency  range,  mainly  in  Arbitrary  Bodies".  Prog.  Aero.  Sciences,  Vol. 

the  sway  coefficients.  This  is  larger  than  the  8. 

corresponding  error  associated  with  equation 

(2.3),  and  is  due  to  the  larger  condition  Kleinaan,  R.  E.  (1982).  "On  the  Mathematical 

number  of  equation  (3.1).  The  effect  of  (3.1)  Theory  of  the  Motion  of  Floating  Bodies  -  An 

in  the  composite  equation  (3.3)  is  here  re-  Update".  DTMSRDC  Report  82/074. 

duced  by  the  selection  of  a  value  for'equal  to 

0.2.  Mei,  C.  C.  (1978).  "Numerical  Methods  in 

Water-Wave  Diffraction  and  Radiation".  Ann. 
‘In  all  cases  tested,  half  of  the  boundary  Rev.  Fluid.  Kech.,  (10),  pp.  393-416. 
of  the  sections  analysed  has  been  approximated 

by  10  straight  segments.  Away  from  the  irregu-  Nestegard,  A.  and  Sclavounos,  P.  D.  (1984).  "A 
lar  frequencies,  the  agreement  between  the  numerical  Solution  of  Two-Dimensional  Deep 

coefficients  obtained  from  equation  (3.3)  for  Water  Wave-Body  Problems".  J.  Ship  Res.,  Vol. 

} -0.2  and  the  numerical  scheme  of  Nestegard  28,  No.  1,  pp.  48-54. 

and  Sclavounos  is  very  good.  Work  is  currently 
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